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Abstract 

This work explores the relation between magnetic islands, long 
range temporal correlations and heat transport. A low order rational 
surface (t = 3/2) was purposely scanned outward through an Electron 
Cyclotron Resonance Heated (ECRH) plasma in the TJ-II stellarator. 
Density turbulence and the poloidal flow velocity (or radial electric 
field) were characterized using a two channel Doppler Reflectometer. 
Simultaneously, the ECRH power was modulated to characterize heat 
transport, using measurements from a 12 channel Electron Cyclotron 
Emission diagnostic. A systematic variation of the poloidal velocity 
was found to be associated with the stationary t = 3/2 magnetic is¬ 
land. Inside from the rational surface, the Hurst coefficient, quantify¬ 
ing the nature of long-range correlations, was found to be significantly 
enhanced. Simultaneously, heat transport was enhanced as well, es¬ 
tablishing a clear link between density fluctuations and anomalous 
heat transport. The variation of the Hurst coefficient was consistent 
with a magnetohydrodynamic turbulence simulation. 


1 



1 Introduction 


In this work, we study the relation between rational surfaces, long range 
temporal correlations and heat transport in the low shear stellarator TJ-II. 
The relationship between rational surfaces and heat transport has been ex¬ 
plored in a number of earlier works, both at TJ-II [HE] and other stellarator 
devices [3]. In tokamaks, low order rational surfaces are often associated 
with Internal Transport Barriers or the reduction of (heat) transport [HE]- 
However, very little is known about the relation between rational surfaces 
and the Hurst coefficient [B], and this work will focus specihcally on this 
issue. The Hurst coefficient is a measure of the degree of long range tem¬ 
poral correlations mm, and hence an important quantiher of anomalous 
transport. 

In the experiments presented here, the position of the rational surfaces is 
slowly scanned outwards by means of the external induction of Ohmic cur¬ 
rent uni, allowing the exploration of various quantities as a function of the 
magnetic flux by diagnostics observing at a hxed position. From Doppler 
reflectometer (DR) measurements, it is found that a rational surface belong¬ 
ing to a low order rational (3/2) is associated with a localized variation of 
the radial electric held. When this rational surface crosses the DR measure¬ 
ment location, the Hurst coefficient obtained from the DR density huctuation 
amplitude is modihed signihcantly. 

In the same discharges, we also characterize heat transport using the 
modulation of centrally deposited Electron Cyclotron Resonance Heating 
(ECRH) combined with radially distributed Electron Cyclotron Emission 
(ECE) temperature measurements. We observe that transport is locally en¬ 
hanced when the DR Hurst coefficient is increased. Eight similar discharges 
are analyzed, showing that the obtained results are systematic and repro¬ 
ducible. 

The reported observations establish a clear link between (a) the modihca- 
tion of density huctuation properties near a rational surface in the gradient 
region, associated with long range memory ehects, and (b) enhanced heat 
transport. 

The interpretation of these observations hnds support from a model m, 
showing that the radial electric held structure and the Hurst coefficient are 
both modihed close to the rational surface. 

This paper is structured as follows. In Section]^ the experimental setup 
is discussed. In Section we present the experimental results. In Section 


2 


we present the modeling results and interpretation. In Section we discuss 
the results, and in Section we draw some conclusions. 
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2 Experimental set-up and techniques used 

In this work, we study discharges heated by Electron Cyclotron Resonant 
Heating (ECRH). In these discharges, the plasma has a relatively low line 
average density of YLf. ~ 0.5 ■ 10^® m“^, so that it is in the electron root state, 
slightly below the critical density of the electron to ion root conhnement 
transition at TJ-II [121 [13]. 

2.1 Modulated ECRH and ECE measurements 

The ECRH system consists of two gyrotrons with a frequency of 53.2 GHz, 
allowing the injection of up to 2 x 300 kW of heating power [Tl|. In the ECRH 
experiments, both ECRH systems were launching into the core of the plasma; 
the half width of the deposition prohle being wecrh — 3 cm [15]. One was 
providing continuous injection of approximately 250 kW, whereas the second 
system was modulated at low power (~ 50 kW), at a modulation frequency 
of 360 Hz and with a duty cycle of 30 %, in order to allow quantifying heat 
transport. 

TJ-H disposes of a 12 channel Electron Cyclotron Emission (ECE) de¬ 
tection system to measure the local electron temperature at up to 12 
different radial positions along the midplane on the high magnetic held side 
of the plasma (at 0 = 315°), covering a signihcant part of the plasma minor 
radius, with a radial resolution of about 1 cm [16]. The local modulation am¬ 
plitude and phase (relative to a central ECE channel) are determined using 
a Fourier analysis of the Te(r, t) data. 

Assuming transport is diffusive it is possible, in principle, to extract dif¬ 
fusive transport parameters (such as the effective diffusivity and convective 
velocity) from these data. However, we refrain from doing so, as the uncer¬ 
tainties introduced by assuming a specihc transport model, as well as the 
ill-posed nature of the problem, do not allow obtaining a clear and unam¬ 
biguous interpretation of the transport processes inside the plasma from such 
an inversion [muH]. Instead, we will represent the amplitude and relative 
phase directly, noting that the gradient of the relative phase is directly linked 
to the speed of the transport [H] l20] . 
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2.2 Doppler reflectometry 

The two-channel Doppler reflectometer of TJ-II [21] is located in a top view¬ 
port at 0 ~ 337°. In Doppler reflectometry, a hnite tilt angle is pnrposely 
introdnced between the incident probing beam and the normal to the reflect¬ 
ing cnt-off layer, and the Bragg back-scattered signal is measnred [22] • The 
amplitnde of the recorded signal is a measure of the intensity of the density 
fluctuations, h. As the plasma rotates in the reflecting plane (a flux sur¬ 
face), the scattered signal experiences a Doppler shift. The size of this shift 
is directly proportional to the rotation velocity of the plasma turbulence 
perpendicular to the magnetic field lines, v±, and therefore to the plasma 
background E x B velocity, provided the latter dominates over the phase 
velocity of density fluctuations [23] . The Doppler reflectometer signals, sam¬ 
pled at 10 MHz, allow determining n and v±_ with high temporal and spatial 
resolution at two radial positions [23] • In the discharges analysed here, chan¬ 
nel 1 is characterized by: p ~ 0.75, ~ 6 cm“^ and channel 2 by: p ~ 0.68, 

fc_i_ ~ 6.7 cm“^. 

2.3 Sweep and tracking of rational surfaces 

In these ECRH experiments, a number of magnetic configurations were used 
with similar (low) magnetic shear, but a slightly different value of rotational 
transform, t. By inducing an Ohmic current in the plasma, the rotational 
transform prohle is modihed, as shown in Fig. [T| The induced Ohmic current 
moves the rational surfaces (such as t = 3/2) outward. Simultaneously, the 
magnetic shear is enhanced, although it remains rather low. 

For the interpretation of experimental data obtained in such discharges, 
a reconstruction of the actual rotational transform profile (and the location 
of the rational surfaces) during the sweep of the Ohmic current is useful. 
This is obtained by iteratively updating equilibria calculated by VMEC [25] . 
using an estimate for the Spitzer resistivity based on experimental data [2S] • 
Details of this calculation will be reported in a separate publication. An 
example result is shown in Fig. with respect to a reference time At = 0 
to be defined in the following. However, this evolution will depend on the 
precise shape of the plasma current prohle, which is not measured directly, 
leading to some uncertainty in the calculation. The interpretation of the 
data reported in the following will therefore not rely on this calculation in 
detail, although it serves as a guide to what one may expect to hnd in the 
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Figure 1: Profiles of the rotational transform for several magnetic configu¬ 
rations and an example of the modihcation produced by an inducted Ohmic 
current. 

experimental situation. 

In this work, we will focus mainly on the t = 3/2 rational and its associ¬ 
ated island. Fig. shows Poincare plots for the poloidal cross sections where 
the ECE diagnostic and the DR diagnostic are located. The t = 3/2 island 
has been generated by adding an appropriate resonant perturbation to the 
vacuum magnetic held corresponding to the magnetic conhguration labelled 
100-36 in Fig. Neither the conhguration, the width of the island, nor its 
precise poloidal location need coincide closely with the actual (dynamically 
varying) experimental situation at any point in time; however, we show the 
hgure as a guide to the interpretation of the measurements reported in the 
following section. 

2.4 The Hurst coefficient 

In this work, we will be using the Hurst coefficient as a measure of long- 
range correlations. This coefficient provides an estimate of the degree to 
which transport deviates from normal dihusive transport. We recall that the 
Hurst coefficient lies in the range 0 < FT < 1. When H = 0.5, one has 
ordinary dihusive transport; when 0 < H < 0.5, one has antipersistence of 
the huctuations, associated with subdihusive transport, and when 0.5 < if < 
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Figure 2: Temporal evolution of the radial location of three rational surfaces 
during the scan of the plasma current (t = 8/5, 3/2, and 7/5). 
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Figure 3: Poincare plots for the vacuum configuration 100_36 (cf. Fig. 0 , in 
0 = constant planes, indicating the positions of the ECE diagnostic (left) 
and the DR diagnostic (right). i?o = 1-5 m. 
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1, persistence of the fluctuations, associated with superdiffusive transport |H1 

The Hurst coefficient is calculated from the density fluctuation amplitude, 
|h|, measured using the DR, using the method described in |2^. We calculate 
the rescaled range graph of R/S using overlapping time windows with a 
length of 10 ms, and evaluate H from the logarithmic slope of R/S{t) over 
the time lag range of 0 < log^g t < 1 (r in ms). The value of H thus obtained 
corresponds to asymptotic times, close to the conhnement time of the plasma. 
Note that the DR sampling rate is 10 MHz, so that 10 ms of data corresponds 
to 10^ data points, providing ample statistics for the determination of the 
Hurst coefficient. Even though the plasma state is varying constantly due to 
the ramped plasma current, these changes are sufficiently slow to consider 
the density fluctuation amplitude measurement to be essentially in steady 
state during the analysis time windows of 10 ms. Modifying the length of 
these time windows by a factor of two hardly affects the results, conhrming 
this statement. Furthermore, MHD mode activity is weak and is found not 
to affect the R/S statistic. 


3 Experiments 


3.1 Doppler reflectometry results 

Figs. 1^ show examples of the analysis of Doppler reflectometry signals. Here, 
v± is compnted from the raw complex signal nsing 64-point bins [2l]. The 
hgnres show v± itself (smoothed to remove noise), showing the mean perpen- 
dicnlar rotation velocity at the observation point, and the Hnrst coefficient, 
H, calcnlated as described in Section 2.4 The 360 Hz modnlation seen in v±_ 


is a conseqnence of the modnlated ECRH power. 

The Hnrst coefficient tends to be slightly above H = 0.5 on average, 
althongh dnring specihc time intervals, the valne is lower or higher, which 
may be related to the passage of rationals by the observation point. In Fig. 

H is particnlarly low for discharge 29782, 1175 < t < 1185 ms, coincident 
with v±_ ~ 0. Similarly, for discharge 29786, H is low for 1205 < t < 1220 
ms, again coincident with n_L ~ 0. The dip in H is somewhat less clear for 
discharge 29787. As will be clarihed in Section]^ the periods of 77 < 0.5 are 
associated with snbdiffnsion occnrring in ‘trapping zones’, characteristic of 
magnetic island 0-points associated with a low order rational snrface. On 
the other hand, in Fig. 77 is high for 1210 < t < 1220 ms in discharge 
29782, and for 1225 < t < 1250 ms in discharges 29786 and 29787. 

Fig. shows the spectrnm of a bolometry channel (ABOL7) for discharge 
29786, also shown in Fig. |2H]- One observes several modes, inclnding a 
mode with a freqnency of ~ 30 kHz, visible in the time interval 1170 < t < 
1220. The mode freqnency is modnlated by the ECRH modnlation freqnency 
(360 Hz). The mean freqnency of the mode sweeps down slowly becanse the 
rational snrface is moving ontward (to parts of the plasma that are rotating 
more slowly) [22] • The apparent disappearance of this mode at f ~ 1220 is 
signihcant, as we will discnss in the following. 

To get a clearer view of the temporal evolntion of the varions qnantities, 
we will nse 8 similar discharges to compnte averages. These discharges are 
listed in Table The starting conhgnrations are different, having vacnnm t 
prohles as shown in Fig. meaning (mainly) that the time of occnrrence of 
the passage of the rational snrface past the Doppler reflectometry observation 
point is different for each discharge. Also, the sweep rates of the OH cnrrent 
vary slightly. To be able to average over snch discharges, we mnst dehne a 
reference time, which we take to be the time of steepest slope of n_L, i.e., the 
times indicated by the vertical dashed lines in Fig. The reference times 
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Figure 4: Smoothed v± signal and Hurst coefficient calculated from \n\ (three 
discharges, DR channel 1). 
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Figure 5: ABOL7 spectrum. The overall slowdown of the frequency is due 
to the rising OH current, slowly moving the rational surfaces outward. The 
modulation of the mode frequency (at 360 Hz) is due to ECRH modulation 
(discharge 29786). 



are different for each channel, as channel 1 (further outward) responds some 
15 ms after channel 2. 

Fig. shows the mean evolution of 77 in a time window centered about 
the reference time. The mean Hurst coefficient of channel 1 is close to the 
value for uncorrelated random noise (diffusion), 77 ~ 0.5, for At < 0, and 
attains a peak of about 77 = 0.65 — 0.7 lasting about 25 ms after the steepest 
descent point of v±. The dip in 77, observed in Fig. is not seen clearly here, 
due to the mentioned variability of the discharge conditions; but the rise in 
77 seen at At > 0 appears to be quite robust to the averaging procedure. 
Channel 2 (further inward) does not show a clear temporal evolution of 77. 

To understand this situation. Fig. shows the mean electron density 
prohles reconstructed from AM reflectometer [HD] and interferometer data, 
using a Bayesian technique described elsewhere [SI]. The density prohles 
evolve only little in the relevant time interval. Clearly, when the rational 
t = 3/2 is located at a position p < 0.65, the density gradient at the rational 
surface is low, whereas it is much higher for p > 0.65. As island growth is 
related to local gradients, we surmise that the corresponding island is initially 
small, but grows when the rational surface reaches the gradient region, p > 
0.65. The rational surface is associated both with a dip in 77 and ~ 0 


11 








Table 1: Discharges and DR reference times 


Discharge number 

Starting configuration 

Reference time 
(Channel 1) 

Reference time 
(Channel 2) 

29771 

100_44_64 

1160 

1150 

29780 

100_46_65 

1170 

1160 

29782 

100_48_65 

1195 

1185 

29783 

100_48_65 

1185 

1175 

29784 

100_50_56 

1210 

1190 

29785 

100_50_56 

1200 

1180 

29786 

100_52_66 

1225 

1200 

29787 

100_52_66 

1225 

1205 



A t (ms) 


Fignre 6: The mean Hurst coefficient (calculated from \h\ for 8 discharges) 
for Doppler reflectometry channels 1 and 2. The red dashed line indicates 
the smoothed trend. 
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Figure 7: Mean electron density profiles for the eight discharges analyzed in 
this work. The profiles are averaged over —100 < At < 100 ms; in this time 
window, temporal variations are only slight. 


(cf. Fig. |4]), the latter being consistent with the stationarity of the island 
in the density gradient region. The stationarity of the island for p > 0.65 
may explain the disappearance of mode oscillations observed by bolometry. 
After passage of the rational, H clearly increases for channel 1 (well inside 
the gradient region at p ~ 0.75), while it remains modest at channel 2 (at 

p ~ 0.68). 


13 















3.2 ECE 


Fig. shows an example of the measured ECE time traces of The 

corresponding power spectra are shown in Fig. The frequency of the hrst 
three harmonics of the modulation frequency (360 Hz) is indicated by vertical 
dashed lines. The amplitude of the hrst three harmonics versus p is shown in 
Fig. [T^ As the modulated ECRH power is deposited centrally, the amplitude 
decays from the centre of the plasma towards the edge. The indicated error 
bars indicate the variation of the amplitude when the analysis is performed 
using successive time windows with a length of 10 ms. The phase of the hrst 
harmonic versus p is shown in Fig. 11, along with a straight line ht. For 


this purpose, a central ECE channel is taken as reference; evidently, for the 
selected reference channel the cross phase is zero. Generally speaking, the 
cross phase increases from the plasma core towards the edge, indicating a 
delayed response to the core modulation. It is understood that the ehect 
of the modulation propagates outward as a consequence of heat transport. 
Note that the phase does not increase monotonically from the centre towards 
the edge; in particular, for p < —0.4 the phase seems to decrease or at least 
remain constant. This suggests that transport is not purely diffusive in this 
region, as pure diffusion would correspond to a monotonically increasing 
phase outside the power deposition region. 

The propagation analysis is quite reproducible for the set of 8 discharges 
cited above. Hence, we use the same reference times as above and calculate 


discharge averaged values versus time and radius. Fig. shows the mean 
modulation amplitude and phase (relative to the central ECE channel). 

Note the signihcant and short-lasting reduction of phase around Af ~ 0 
ms (IpI ~ 0.65). This time and place coincides rather well with the passage 
of the rational t = 3/2, cf. Figj^ and its entry into the density gradient zone. 
To interpret this behavior, recall that the measured modulation amplitude 
and phase can be related to the effective heat diffusion coefficient via [19]: 


Xe = 


OJq 


/A + l/2r + n'/2n) 


( 1 ) 


Thus, a drop in the radial derivative of the phase, </', implies an increase in 
local X- Fig- [12] shows that the phase is nearly constant over a wide area. 


namely 0.3 < p < 0.7, for a time period lasting several tens of ms around 
At = 0. This observed reduction of modulation phase gradient therefore 
indicates a spatially and temporally localized enhancement of heat transport. 
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Figure 8: ECE data Tf,{p,t) obtained in an ECRH modulation experiment 
(discharge 29786, ECE channels 1,3,6,8,11). 
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Figure 10: Amplitude, A, of the first three harmonics versus radius, calcu¬ 
lated from the Te(p, t) data shown in Fig. 
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Figure 11: Cross phase, 0, of the hrst harmonic versus radius, calculated 
from the Te(p, t) data shown in Fig. 
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Figure 12: Mean ECE modulation results (top: A, bottom: 0). 
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The slight temporal variation of A is mainly due to a small temporal variation 
of n^, leading to a corresponding small variation in absorbed ECRH power. 

Fig. shows the time evolution of the measured electron temperature, 
(Te(p, f)), averaged over the 8 mentioned discharges, and its temporal varia¬ 
tion, A(Te(p, t)) = (Te(p, f)) — (Te(p, f)), where the overline indicates a time 
average taken over the time window —50 < Af < 0 ms. By dehnition, this 
quantity is near zero for —50 < At < 0 ms, but a significant temperature 
drop (|ATe| < 50 eV) is visible for At > 0, starting at around p = 0.35 and 
propagating both inward and outward, suggesting a release of heat from the 
core towards the edge. 
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Figure 13: {Te{p,t)), averaged over 8 similar discharges, and A(Te(p,t)). 
Units: keV. 
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3.3 ECE modulation amplitude variations 


Fig. 14 shows the signals of the outermost 4 ECE channels (top). The mod¬ 


ulation signal ‘disappears’ in ECE channel 2, for 1165 < t < 1188 ms, the 
time interval indicated by the vertical dashed lines. The length of this time 
interval coincides with the half period of the velocity variation observed by 
DR (not shown for this discharge, but similar in shape to Fig. discharge 
29782), which in turn is related to the radial extension of the magnetic is¬ 
land [32l [33] . Thus, a possible explanation for this phenomenon is that the 
modulated component of is weak inside the island 0-point region as it 
sweeps past the observation point. Even so, the modulation of the centrally 
deposited ECRH power does reach adjacent ECE channels, so this effect oc¬ 
curs only at specihc radii. No important prohle flattening is observed, within 
the resolution of the diagnostic |M|. Thus, the reason for the disappearance 
of the modulated component could be related to the magnetic topology, the 
island being isolated from the main plasma, and/or the relatively strong 
sheared flow associated with the rational, in combination with a signihcant 
local gradient. 

Again, to show how systematic this reduction of the modulated com¬ 
ponent of Te is, we calculate an average over the 8 mentioned discharges. 


For this purpose. Fig. [^displays the mean ratio R, a quantiher of the de¬ 
gree of isolation from the temperature modulation, versus time and radius. 
This ratio R = R 1 /R 2 is dehned as follows. Ri is dehned as the root mean 
square (RMS) amplitude of the measured Te{p,t) data, after high-pass hlter- 
ing with a cutoff frequency of 2 kHz. R 2 is dehned as the RMS amplitude of 
the measured Te{p,t) data, after band-pass hltering in the frequency range 
0.1 < / < 2 kHz. Thus, Ri is a measure of the noise amplitude while R 2 is a 
measure of the modulation amplitude (the modulation frequency being 0.36 
kHz). 

As argued above, the island 0-point region is topologically separate from 
the main plasma and/or characterized by a signihcant sheared how (at least, 
in the present case), which would reduce the modulated component of 
in this region. Thus, the regions in Fig. 15 with high values of R (i.e., low 


modulation amplitude) would correspond to island 0-point regions and/or 
sheared how regions. Such regions with high R can of course only be observed 
clearly (at the hxed ECE measurement positions) if the island is large and 
stationary, and/or the sheared how is strong. Apparently, this condition 
is met when the rational 3/2 surface (moving outward) reaches |p| > 0.65, 
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Figure 14: Measured ECE temperature (discharge 29780). Shown are ECE 
channels ECEl through ECE4, at p = —0.75, —0.64, —0.53, —0.44, respec¬ 
tively. Units: keV. 
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Figure 15: Ratio R of high frequency noise versus modulation amplitude, 
measured by the ECE system (see text), plotted versus time and radius 
(averaged over 8 discharges). 
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which is the density gradient region. At earlier times, when the rational is 
farther inward, the island may be smaller and/or rotating, snch that R is 
smaller. The zone of relatively high R seen at |p| ~ 0.4 is possibly associated 
with another rational (t = 7/5); cf. Fig. In view of the above argnment, 
this 7/5 mode mnst also be stationary to be visible; yet it is not located in 
the gradient region. Therefore, we specnlate that its apparent stationarity 
is a conseqnence of mode locking with the stationary 3/2 mode. Note that 
the total magnetic field, \B\, varies by less than 0.5% dne to the varying 
plasma cnrrent, so that the radial locations of the ECE channels, inversely 
proportional to \B\, may be considered fixed. 


25 


4 Modeling 

Topological analysis of the flow structures of plasma turbulence [35] has 
shown that typically, low order rational surfaces are associated with small 
transport barriers. These are created by zonal flows, induced by flow eddies 
associated with the rational surfaces. These small transport barriers tend 
to trap particles, leading to local subdiffusive behavior, as reported in [36] . 
Therefore, a reduction of the Hurst coefficient below 0.5 is expected in the 
barrier regions. On the other hand, experimental studies have shown that 
H is systematically large in the plasma edge region [S] and varies with the 
global gradients EH. 

To study the correlation between the zonal flows associated with low 
order rationals and the Hurst coefficient, we have performed a simulation 
of the magnetohydrodynamic turbulence induced by resistive interchange 
modes in a periodic cylinder. Details of this type of modeling can be found 
in Ref. m- We have not attempted to match the experimental situation 
closely, due to the difficulty of reproducing a gradual modihcation of the 
rotational transform in time, induced by an increasing plasma current. Here, 
we limit ourselves to studying the possible relation between low order rational 
surfaces, zonal flow generation and the Hurst coefficient. For this purpose, 
we have taken one of the standard rotational transform prohles of TJ-H, 
the case labelled 100_46 in Fig. [T| using a pressure prohle that is close to 
the experimental one. Because the rotational transform has low shear, the 
low order rational surfaces produce magnetic islands of a signihcant width, 
facilitating the detection of regions with potential subdiffusion. 

The main parameters used in these calculations are: /3q = 0.001, inverse 
aspect ratio e = a/R q = 0.15 and Lundquist number S = tr/ta = 10®. 
For these parameters, this conhguration is unstable to resistive interchange 
modes. Here, /3o is the ratio of the plasma pressure and the magnetic pressure 
at the plasma axis, (3q = 2/iop(0)/R^, tr is the resistive time at the magnetic 
axis, Tr = /r]{Q), where ri{0) is the resistivity at the magnetic axis and 

ta is the Alfven time, ta = RQ^HomiTii/B(^. 

We have calculated the Hurst coefficient at 100 equally spaced radial 
locations over a period of time of Atr, long enough to obtain asymptotic 
behavior. We also calculated the angular and time averaged poloidal flow. 
The Hurst coefficient and the averaged poloidal flow are shown in Fig. 

At each of the main rational surfaces, we can see that the averaged 
poloidal flow has a minimum, while the flow has strong shear in the region 
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Figure 16: Modeling results: mean poloidal flow and Hurst coefficient near 
various rationals. 

of the magnetic island associated with each resonant surface. It is also clear 
that the H coefficient has a minimum in the same region, with values below 
0.5. The H coefficient increases above 0.5 as we move out of the region of the 
resonant surface. This structure is fairly systematic at all the low resonant 
surfaces present in the plasma and they are consistent with the measurements 
shown in Fig. The experiment also shows a clear enhancement of H inside 
from the rational surface, cf. Fig. 

When comparing with the experiment, on should remember that trans¬ 
port consists of both turbulent and collisional transport, the latter being 
induced by collisions |3H]. Collisional transport has a diffusive nature and 
tends to dominate subdiffusive processes, making it difficult, in general, to 


27 
















r 


Figure 17: Modeling results: instantaneous poloidal flows at two time points, 
and mean poloidal flow. 


observe significant regions with H < 0.5 in the experiment. Nevertheless, 
times with low values of H are visible in Fig. and H is systematically 
moderate (of the order of 0.5 or slightly below) in a zone associated with the 
t = 3/2 island, cf. Fig. 

(ye), is calculated by averaging the 


The averaged flow shown in Fig. 16 


poloidal flow over the poloidal and toroidal angles and over a period of time 
such that the dynamical calculation and the corresponding zonal flow are in 
steady state. Note that this mean flow is different from the perpendicular 
flow measured in the experiment by DR, which is instantaneous and local, 
and hence may be very sensitive to the precise poloidal and toroidal location 
of the measurement (i.e., the phase of the stationary island). If we compare 
the averaged flow near a singular surface to the instantaneous poloidal flow, 
cf. Fig. one observes that the latter is not symmetric with respect to 


the singular surface, and that the asymmetry may vary. This would explain 
the asymmetric appearance of the poloidal flow with respect to the rational 
surface observed in the experiment. 

The H coefficient is determined by at least two effects. One is the gen- 
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Figure 18: Modeling results: absolute value of the gradient and Hurst coef¬ 
ficient. 


eration of long range correlations by avalanches, and the other is their sup¬ 
pression due to decorrelation by sheared flow. Avalanches are non-local phe¬ 
nomena, but in the plasma core region avalanches are relatively rare and 
they tend to correlate with the gradient, as they originate there. The edge, 
however, is crossed by many avalanches initiated in the entire plasma, so no 
relation is expected between the local gradient and the H exponent. At the 
plasma edge, the H exponent tends to be above 0.5 [8]. In Fig. 18 we show 
the gradient of the pressure, averaged over the flux surface and over time, 
as a function of the radius, together with the H coefficient, to illustrate this 
situation. At the singular surface, the density gradient is low, and due to 
the presence of the islands and the associated average sheared flow, the H 
coefficient is reduced below 0.5. Immediately outside the islands, where the 
sheared flow goes to zero, the density gradient is high, correlated with an 
increase of H above 0.5, at least in the plasma core region. In the plasma 
edge region, H tends to rise above 0.5 globally, while it is reduced below 0.5 
only at the singular surfaces, mainly t = 5/3 in this case, without any clear 
correlation with the gradient. 
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5 Discussion 


Doppler reflectometry (DR) has detected a structure in v± (equivalent to Er) 
that is systematically associated with the t = ?)/2 rational surface, which is 
purposely and slowly scanned ontward in the experiments discussed here. 
We take the time of maximum temporal variation of v±_ as a reference time 
(At = 0); this time is closely associated with, though not identical to, the 
time of passage of the rational snrface itself |33]. When the rational surface 
crosses the DR location of channel 1, the Hnrst coefficient, H, is rednced (cf. 
Fig. 0. This is consistent with the systematic reduction of H at rational 
snrfaces observed in the model discussed in Section IH A rednction of H 
is associated with subdiffusive transport, implying that possibly, the island 
constitntes a (minor) transport barrier [391110], which in tnrn may be related 


to the relatively good heat confinement observed in Fig. 13 for At < 0. 

If the location of the singular surface is associated with the minimum of 
the Hurst coefficient. Fig. [^indicate that the passage of the singular surface 
occurs slightly before (i.e., it’s outward) from the point of maximum slope of 

U_L. 

Some MHD mode oscillations (~ 30 kHz) are observed (bolometry, DR) 
for At < 0, which are rednced or disappear afterwards. As these mode 
oscillations gradually reduce their mean freqnency, it is clear that they are 
associated with rational snrfaces that are moving ontward to areas of the 
plasma with smaller mean rotation velocity [29]. The proposed explanation 
for the disappearance is that the t = 3/2 island is small for At < 0, bnt large 
and stationary for At > 0. This is consistent with the observation u_i_ ~ 0 at 
the rational surface in the density gradient region, cf. Fig. [^ 

This interpretation, which relies on the existence of a large island in the 
gradient region, hnds snpport from the observation of specihc radial zones 
where the ECRH modulation penetrates only weakly (Fig. 15), which are 
interpreted as topologically isolated island 0-point regions [iQ] and/or zones 
with considerable velocity shear. 

For 0 < At < 25 ms, the Hurst coefficient is large {H{p ~ 0.75) ~ 
0.65 — 0.7), althongh somewhat farther inward, H stays moderate {H{p ~ 
0.68) ~ 0.5). A representation of the evolntion of the mean ECE temperature 
(Tg) and its variation A(Te) (Fig. 13) snggests heat loss occurring for At > 
0, starting in the core and expanding ontward, as visnalized in the plot 
of A(Te) as a dark bine area. Simultaneously, the Tg modulation phase 
delay gradient (d^/dr) is reduced in the region 0.3 < p < 0.7. We note 
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that the (incremental) heat diffusivity is proportional to the inverse of the 
gradient of the phase m, so that the decrease of d(j)/dr suggests enhanced 
heat transport and stiffness in this region (inward from the island). This 
observation establishes a clear link between (memory effects associated with) 
density fluctuations and heat transport. 
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6 Conclusions 


In these experiments at TJ-II, an externally driven OH current was used to 
modify the rotational transform prohle, t(p), and slowly move the rational 
surfaces outward through the plasma. As the main rational t = 3/2 moved 
past the DR observation points, the perpendicular propagation velocity of 
the turbulence, v± - equivalent to the radial electric held, Ey - was found 
to experience a signihcant ‘wiggle’. Taking the time of strongest temporal 
variation of v± as a reference, 8 similar discharges were analyzed. 

Small rotating MHD modes, were detected by bolometry, with decreasing 
frequency, consistent with the radial outward motion of the rational surface 
locations. Once the t = 3/2 rational was well inside the density gradient 
region, the corresponding island grew and was found to be stationary. No 
clear battening of prohles was observed, which could be consistent with ex¬ 
pectations for small magnetic shear [32]. However, it was found that the 
modulated heat deposited in the core plasma only penetrated weakly in the 
zone associated with the stationary island. 

The Hurst coefficient was found to be low [H < 0.5) at the rational 
surface (Fig. and systematically high [H ~ 0.65 — 0.7) inward from the 
t = 3/2 rational surface (Fig. [^, as long as the rational surface was located 
well inside the density gradient region. The analysis of ECE data revealed 
that radial heat transport was enhanced in the region inward from the large 
island. 

Simulations of magnetohydrodynamic turbulence induced by resistive in¬ 
terchange modes reproduced some of these observations qualitatively, namely: 
(1) the variation of poloidal how prohles associated with rational surfaces, 
and (2) the variation of the Hurst coefficient associated with rational surfaces 
and local pressure gradients. 

Thus, this work establishes a clear correlation between density huctu- 
ation long range temporal correlation properties (quantified via the Hurst 
coefficient) and enhanced heat transport in the presence of a large island, 
possibly providing the hrst direct observation of the turbulent mechanism 
behind anomalous heat transport. 
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